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In the light of recent experimental results for the direct detection of dark matter, we 
analyze in the framework of SUGRA the value of the neutralino-nucleon cross section. We 
study how this value is modified when the usual assumptions of universal soft terms and 
GUT scale are relaxed. In particular we consider scenarios with non-universal scalar and 
gaugino masses and scenarios with intermediate unification scale. We also study superstring 
constructions with D-branes, where a combination of the above two scenarios arises naturally. 
In the analysis we take into account the most recent experimental constraints, such as the 
lower bound on the Higgs mass, the b — > sj branching ratio, and the muon g — 2. 

1 Introduction 

One of the most interesting candidates for dark matter is a long-lived or stable weakly- 
interacting massive particle (WIMP). WIMPs can remain from the earliest moments 
of the Universe in sufficient number to account for a significant fraction of relic density. 
These particles would form not only a background density in the Universe, but also 
would cluster gravitationally with ordinary stars in the galactic halos. This raises the 
hope of detecting relic WIMPs directly, by observing their elastic scattering on target 
nuclei through nuclear recoils [|TJ. In fact, one of the current experiments, the DAMA 
collaboration, reported recently data favoring the existence of a WIMP signal. It 
was claimed that (at 4a C.L.) the preferred range of the WIMP-nucleon cross section 
is 10~ 6 -10~ 5 pb for a WIMP mass between 30 and 200 GeV. Unlike this spectacular 
result, other collaborations, CDMS §, EDELWEISS §, and IGEX §, claim to have 
excluded regions of the DAMA parameter space. 

In any case, due to these and other projected experiments, it seems very plausible 
that the dark matter will be found in the near future. In this situation, and assuming 
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that the dark matter is a WIMP, it is natural to wonder how big the cross section 
for its direct detection can be. The answer to this question depends on the particular 
WIMP considered. The leading candidate in this class is the lightest neutralino [J], 
a particle predicted by the supersymmetric (SUSY) extension of the standard model 
(SM). 

In the simplest SUSY extension, the so-called minimal supersymmetric standard 
model (MSSM), there are four neutralinos, x° {i = 1, 2, 3, 4), since they are the physical 
superpositions of the bino (B°), of the neutral wino (W®), and of the neutral Higgsinos 
(H®, Hj). In most of the parameter space of the MSSM, the lightest neutralino, Xii 
is the lightest supersymmetric particle (LSP). When R-parity is imposed this implies 
that Xi is absolutely stable, and therefore a dark matter candidate. 

In this paper we will analyze this SUSY scenario in the framework of supergravity 
(SUGRA), and in particular several constructions proposed recently in the literature 
(for a review see ref. ||), where the Xi-nucleon cross section can be enhanced, and might 
be of the order of 10~ 6 pb, i.e., where current dark matter detectors are sensitive. First, 
in Section 2, we will briefly introduce the SUGRA framework. In Section 3 we will 
describe the most recent experimental constraints which can affect the computation 
of the cross section. Then, in the rest of the sections, we will re-evaluate previous 
computations, taking into account these constraints. In particular, in Section 4 we 
will review the value of the cross section when universal soft SUSY-breaking terms are 
assumed, and how this value is modified when this assumption is relaxed f7|, || |, || |10| . 
In Section 5 we will consider the case of an intermediate unification scale [O, IT3. 



Finally, in Section 6, we will study superstring scenarios with D-branes [13 



2 SUGRA framework 

Working in the framework of SUGRA one makes several assumptions. In particular, 
the soft parameters, i.e., gaugino masses, scalar masses, and trilinear couplings, are 
generated once SUSY is broken through gravitational interactions. They are denoted 
at the grand unification scale, Mqut ~ 2 x 10 16 GeV, by M a , m a , and A a ^ respectively. 
Likewise, radiative electroweak symmetry breaking is imposed, i.e., the Higgsino mass 
parameter fi is determined by the minimization of the Higgs effective potential. This 
implies 

mir — mi, tan 2 3 1 
M = \ - n M Z , (i) 



where tan (3 = (H®) / (H%) is the ratio of Higgs vacuum expectation values. 

The soft SUSY-breaking terms will have in general a non-universal structure in 
the framework of SUGRA |TJ]]. For the case of the observable scalar masses, this is 
due to the non-universal couplings in the Kahler potential between the hidden sector 
fields breaking SUSY and the observable sector fields. For the case of the gaugino 
masses, this is due to the non-universality of the gauge kinetic functions associated 
to the different gauge groups. Explicit string constructions, whose low-energy limit is 
SUGRA, exhibit these properties . 
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With these assumptions, the SUGRA framework still allows a large number of free 
parameters. In particular, it contains more than 100 parameters. In order to have 
predictive power one usually assumes in addition that the above soft parameters are 
real and universal at Mqut- This is the so-called minimal supergravity (mSUGRA) 
scenario, where there are only four free parameters: m, M, A, and tan f3. In addition, 
the sign of fi remains also undetermined. It is worth noticing that explicit string 
constructions with these characteristics can also be found [|i~4|. 



3 Experimental constraints 

Here we list the most recent experimental results which might be relevant when com- 
puting the x^-nucleon cross section in the context of SUGRA. In particular, they may 
produce important constraints in the parameter space, forbidding for example points 
with a large cross section. 

• Higgs mass 

Whereas in the context of the SM, the negative direct search for the Higgs at 
the LEP2 collider implies a precise lower bound on its mass of 114.1 GeV, the 
situation in SUSY scenarios is more involved. In particular, in the framework 
of mSUGRA, one obtains [[DJ for the lightest CP-even Higgs mass rrih > 114.1 
GeV when tan/3 < 50, and > 91 GeV when tan/3 is larger. Let us remark 
that tan/3 is constrained to be smaller than 60 since otherwise the condition of 
radiative electroweak symmetry breaking cannot be fulfilled. 

On the other hand, when the mSUGRA framework is relaxed the above SUSY 
bounds must be modified. In particular, for the benchmark scenarios studied in 
ref. |nj one obtains > 114.1 GeV when tan/3 < 8, and m,h > 91 GeV when 



tan (3 is larger. 

Let us finally remark that in our computation below we evaluate using the 
program FeynHiggsFast, a simplified version of the updated program FeynHiggs 



17J which contains the complete one-loop and dominant two-loop corrections. 
The value of obtained with FeynHiggsFast is approximately 1 GeV below the 
one obtained using FeynHiggs. In addition, we should also keep in mind that the 
value of irth obtained with FeynHiggs has an uncertainty of about 3 GeV, due 
e.g. to higher-order corrections. 

Top mass 

Needless to say we use as input for the top mass throughout this paper the 
central experimental value mt(pole) = 175 GeV. However, let us remark that a 
modification in this mass by ±1 GeV implies, basically, a modification also of ±1 
GeV in m^. 

SUSY spectrum 
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We impose the present experimental lower bounds on SUSY masses coming from 
LEP and Tevatron. In particular, using the low-energy relation from mSUGRA, 



Mi = | tan 9wM2, one obtains for the lightest chargino mass the bound |18 



m~± > 103 GeV. Likewise, one is also able to obtain the following bounds for 



sleptons masses [0 : mg > 99 GeV, > 96 GeV, m f > 87 GeV. Finally, we use 
the following bounds on the masses of sneutrino, the stop, the rest of squarks, 
and gluinos: rxiy > 50 GeV, m t - > 95 GeV, m- q > 150 GeV, m- g > 190 GeV. 

b — > S7 



The measurements of B — > Xg-y decays at CLEO [^0J and BELLE pi] , lead to 
bounds on the branching ratio b —>■ S7. In particular we impose in our computa- 
tion 2 x 10~ 4 < BR(b —>■ S7) < 4.1 x 10~ 4 , where the evaluation is carried out 



at the NLO accuracy in the SM [2~2~|, while for the SUSY contributions only the 



1-loop diagrams have been taken into account |2 

The new measurement of the anomalous magnetic moment of the muon, a M = 
(g^ — 2)/2, in the E821 experiment at the BNL implied apparently a 2.6a 
deviation from the SM predictions. In particular, taking a 2a range around 
the E821 central value one would have 11 x 10~ 10 < a M (£821) - a^(SM) < 
75 x 10~ 10 . However, recent theoretical computations |25| have shown that a 
significant part of this discrepancy was due to the evaluation of the hadronic 
light-by-light contribution p6| . As a consequence, assuming that the possible new 
physics is due to SUGRA, we impose in our computation the current constraint 
-7 x 10~ 10 < a^SUGRA) < 57 x lO" 10 at the 2a level. 

LSP 

The LSP must be an electrically neutral (also with no strong interactions) parti- 
cle, since otherwise it would bind to nuclei and would be excluded as a candidate 
for dark matter from unsuccessful searches for exotic heavy isotopes |27] . As men- 
tioned in the Introduction, the Xi is the LSP in most of the parameter space of 
SUGRA. However, in some regions the stau can be lighter than the Xv Following 
the above argument, we discard these regions. 

Other cosmological constraints, as e.g. the observational bounds on the relic Xi 
density, 0.1 < Q^oh 2 < 0.3, will not be applied. The computation of the relic Xi 
density depends on assumptions about the evolution of the early Universe, and 



therefore different cosmological scenarios give rise to different results [28 . 



4 SUGRA predictions for the neutralino-nucleon 
cross section 



The cross section for the elastic scattering of Xi 011 nucleons has been examined ex- 
haustively in the literature (for a recent re-evaluation see ref. [^|). Let us recall that 
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the relevant effective Lagrangian describing the elastic scattering of Xi on protons and 
neutrons has a spin-independent (scalar) interaction and a spin-dependent interaction. 
However, the contribution of the scalar one is generically larger and therefore we will 
concentrate on it. This scalar interaction includes contributions from squark exchange 
and neutral Higgs exchange. Let us also remark that the scalar cross sections for both, 
protons and neutrons, are basically equal. 

In what follows we will review the possible value of the cross section in the frame- 
work of SUGRA. First we will consider the mSUGRA scenario, where the soft terms are 
assumed to be universal, and we will later relax this condition allowing non-universal 
soft masses. 



4.1 mSUGRA scenario 

In the mSUGRA scenario Xi is mainly bino. This result can be qualitatively understood 
from the well known evolution of m 2 Hu with the scale. It becomes large and negative 
at low scales. Then given by eq. (j|D becomes also large, and in particular much 
larger than M\ and Mi. Thus, the lightest neutralino will be mainly gaugino, and 
in particular bino, since, as mentioned in the previous section, at low energy M\ = 
§ tan 2 6 W M 2 w 0.5M 2 . 

When Xi is basically bino scattering channels through Higgs (h, H) exchange are 
highly suppressed, and as a consequence, the predicted scalar x°-proton cross section 
is well below the accessible experimental regions for low and moderate values of tan f3 
29 1 . Although the cross section increases entering in the DAMA region || |3l| when 



the value of tan/3 increases [|7|, |8|], the experimental constraints discussed in Section 3 
exclude this possibility [30, 1C]. We show this fact in Fig. |I|. There the x^-proton cross 



section function of the neutralino mass m^o for tan /3 = 20 and 35 is plotted. 

Concerning the parameter space of the figure, we use the following values. For 
the soft scalar mass 0< m < 600 GeV, where the curve associated to m = 600 GeV 
corresponds to the minimum values of cr^o_ p in the figure. Note that in order to avoid 
the stau being the LSP not very small values of m are in fact allowed. The cross 
section is not very sensitive to the specific values of A in a wide range. In particular, 
we have checked that this is so for | A/M | < 4. In the figure we fix A = M, 
corresponding to the maximum value of cr^o-p- O n the other hand, the gaugino mass 
M is essentially fixed for a given tan/5 and rn^o. We consider only positive values 
of M, since negative ones correspond to smaller cross sections. Finally, we choose to 
plot only the case \i > since for negative values of \i the cross section is smaller. 
In addition, constraints coming from the b — > 37 and — 2 processes highly reduce 
the n < parameter space. Let us remark that we are using the conventions for the 
SUGRA parameters entering in the Lagrangian explained in ref. 0. 

Concerning the experimental constraints, the lower bound in the figure m^o > 50 
GeV is obtained imposing the experimental bound on the lightest chargino mass m-± > 
90 GeV. However, as mentioned in Section 3, the most recent bound in the context of 



mSUGRA is now 181 m-± > 103 GeV, and we show it with a dotted line. We also 



show in Fig. 1 several values of the Higgs mass with dashed lines. In addition the 
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Figure 1: Scatter plot of the neutralino-proton cross section <J^o_ p as a function of 
the neutralino mass m^o in the mSUGRA scenario, for tan/3 = 20 and 35. The dotted 
line corresponds to the LEP lower bound m~± = 103 GeV. Dashed lines correspond 
to different values of the Higgs mass, 110 and 114 GeV. The region to the left of the 
dot-dashed curve is excluded by b — > 57, and the region above the solid line is excluded 
by g M — 2. DAMA and CDMS current experimental limits and projected GENIUS 
limits are shown. 



region to the left of the dot-dashed curve is excluded by b — > 37, and the region above 
the solid line is excluded by — 2. 

As we can see, the present experimental lower limit for the Higgs mass in mSUGRA 
when tan/3 < 50, rrih > 114.1 GeV, is the most important constraint, implying that 
only the solid region where c^o- p ~ 10~ 8 pb survives. Obviously, in this mSUGRA 
case, more sensitive detectors producing further data are needed. Fortunately, many 
dark matter detectors are being projected. Particularly interesting is the GENIUS 
detector [^|, where values of the cross section as low as 10 -9 pb will be accessible, as 
shown in the figure. 



Very large values of tan /3, like tan (3 ~ 50, have also been considered |33|, [!(], |T0 
Although it was found, as expected, that the cross section is enhanced, the well known 
experimental limits coming from b — > 37 for large tan /3, lead to severe constraints on 
the parameter space. In particular, these constraints imply cr^o_ p < 10~ 7 pb. 



4.2 Scenario with non-universal soft terms 

It was shown recently that non-universality allows to increase the neutralino-proton 
cross section. The key point consists of reducing the value of Following the discus- 
sion in the previous subsection, the smaller is, the larger the Higgsino components 
of the lightest neutralino become. Eventually, will be of the order of Mi, M2 and 
Xi will be a mixed Higgsino-gaugino state. Indeed scattering channels through Higgs 
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Figure 2: Scatter plot of the neutralino-proton cross section cr^o_ p as a function of the 
neutralino mass m^o in the SUGRA scenario with non-universal soft scalar masses, for 
tan /3 — 20 and 35. Only points of the parameter space fulfilling b — > 37 and g M — 2 are 
plotted. Big (red) dots correspond to points with rrih > 114 GeV. Small (blue) dots 
correspond to points with 91 < rrih < 114 GeV. The universal case (green region) is 
shown for comparison. DAMA and CDMS current experimental limits and projected 
GENIUS limits are also shown. 



exchange are now important and their contributions to the cross section will increase 
it. This can be carried out with non- universal scalar masses |7], |8|, [K| and/or gaugino 
masses || 

We can have a qualitative understanding from the following. When m 2 H u { m2 H d ) a t 
Mqut increases(decreases), its negative (positive) contribution at low energy in eq. ([!]) 
is less important. Likewise, when and mf lR at M GUT decrease, due to their con- 
tribution proportional to the top Yukawa coupling in the RGE of m 2 Hu the negative 
contribution of the latter to /i 2 is again less important. 

This effect can be seen in Fig. |2], where a scatter plot of function of m^o 

for a scanning of the soft scalar masses is shown. This scanning is explained in detail 
e.g. in ref. ||. As in Fig. [I], we are taking A = M. For comparison we superimpose 
also the region (green area) obtained in Fig. [I] with universality. We see that non- 
universal scalar masses help in increasing the value of cr^o_ p . In the figure we only plot 
the points of the parameter space fulfilling b — > 37 and — 2. In particular, the big 
(red) dots correspond to points with rrih > 114 GeV. As a conservative approach, we 
also plot with small (blue) dots those points with 91 < < 114 GeV (as discussed 
in Section 3, this possibility may arise when the mSUGRA scenario is relaxed). For 
tan/3 = 20, both kind of points enter in the DAMA region. In particular, those with 
rrih < 114 GeV have an upper bound for the cross section w 10~ 5 pb, whereas for those 
with rrih > 114 GeV the upper bound is ~ 10~ 6 pb. For tan/3 = 35, the area of points 
entering in the DAMA region increases a lot, and now, even points with rrih > 114 
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Figure 3: The same as in Fig. [| but using non-universal soft gaugino masses. 



GeV can have c"^o_ p ~ 1CU 5 pb. Let us finally remark that the non-universality in the 
Higgs sector gives the most important effect, and the non-universality in the sfermion 
sector only increases slightly the cross section. 

On the other hand, concerning gaugino masses, it is worth noticing that M 3 appears 
in the RGEs of squark masses, so e.g. their contribution proportional to the top Yukawa 
coupling in the RGE of rn 2 Hu will do this less negative if M3 is small, and therefore /x 2 
will become smaller in this case. Taking into account this effect, we show in Fig. 0, for 
tan/3 = 20 and 35, a scatter plot of o^'-p as a function of m^o for the scanning of the 
parameters M 3 , M 2 , M\ explained in ref. ||. Although points with a large cross section 
entering in the DAMA region do exist ||, they are, at the end of the day, forbidden 
when b — > 57 and — 2 constraints are imposed. 

Clearly, in this sense, SUGRA scenarios with non-universal scalars are favored with 
respect to scenarios with non-universal gauginos. 



5 Scenario with intermediate unification scale 

The analyses of cr^o_ p in SUGRA, described in the previous section, were performed 
assuming the unification scale Mqjjt ~ 10 16 GeV. However, there are several interesting 
arguments in favor of SUGRA scenarios with scales Mi ~ io 10-14 GeV. These can be 
found e.g. in ref. and references therein. Thus to use the value of the initial scale, say 
Mi, as a free parameter for the running of the soft terms is particularly interesting. In 



this sense, it was recently pointed out [TT| that cr^o_ p is very sensitive to the variation of 
the initial scale for the running of the soft terms. For instance, by taking M/ = 10 10-12 
GeV rather than Mqut, regions in the parameter space of mSUGRA can be found 
where cr^o_ p is two orders of magnitude larger than for M GUT . 

The fact that smaller initial scales imply a larger cr^o_ p can be understood from the 
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Figure 4: The same as in Fig. [1] but including the case with intermediate scale, 
Mi = 10 11 GeV. Now the Higgs masses 105 and 110 GeV are shown. 



variation in the value of \i with Mj. One observes that, for tan/3 fixed, the smaller the 
initial scale for the running is, the smaller the numerator in the first piece of eq.(|l]) 
becomes. This can be understood qualitatively from the well known evolution of m 2 Hd 
and tn 2 H with the scale. Clearly, the smaller the initial scale is, the shorter the running 
becomes. As a consequence, also the less important the positive(negative) contribution 
Trig (rriffj) to /i 2 in eq.([Lj) becomes. Thus \fi\ decreases, and therefore, as discussed in 
the previous section, cr^o_ p increases. This is shown in Fig. |], where the result for the 
scale Mj = Mqut is compared with the result for the intermediate scale Mi = 10 11 
GeV. 

Clearly, in the latter scenario the Higgs mass decreases and all points would be 
excluded if we impose the mSUGRA bound rrih > 114 GeV. On the other hand, 
taking into account the conservative approach discussed in Section 3, since all points 
in the figure have > 100 GeV we will not discard them. Then, the most important 
constraint for tan (3 = 20 is due to — 2, excluding a large area of points in the DAMA 



region. We can understand this from the results of ref. There, a M versus m^o is 
plotted, and one can observe that for m^o < 150 GeV and m small too large values 
of a M , beyond the experimental upper bound, are obtained. Note that the area left in 
DAMA by this constraint and the mSUGRA one m-± > 103 GeV may be increased if 
we relaxed the latter bound to e.g. 90 GeV. We see in Fig. |] that the range 75 GeV 
< m^o < 125 GeV is now consistent with DAMA limits. However, for tan/3 = 35 all 
points in the DAMA region are excluded because of b — > sj. 

Given the above situation concerning the enhancement of the neutralino-proton 
cross section through non-universal scalars and intermediate scales, it is worth analyz- 
ing the combination of both possibilities. This is shown in Fig. [5| Although no points 
with rrih > 114 GeV enter in the DAMA region, many points with 91 < < 114 
GeV do enter and are allowed by all the other constraints. For tan /3 = 35 the allowed 
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Figure 5: Scatter plot of the neutralino-proton cross section cr^o_ p as a function of the 
neutralino mass m^o in the SUGRA scenario with the intermediate scale Mr = 10 11 
GeV and non-universal soft scalar masses, for tan/3 = 20 and 35. Only points of 
the parameter space fulfilling b — > 57 and — 2 are plotted. (Blue) dots correspond 
to points with 91 < rrih < 114 GeV. The universal case (green region) is shown for 
comparison. DAMA and CDMS current experimental limits and projected GENIUS 
limits are also shown. 



area is extremely large and we can even have points entering in the region excluded by 
CDMS. 



6 Superstring scenario with D-branes 

In the previous section the analyses were performed assuming intermediate unification 
scales. In fact, this situation can be inspired by superstring theories, since it was 
recently realized that the string scale may be anywhere between the weak and the 
Plank scale. For example, embedding the SM inside D3-branes in type I strings, the 
string scale is given by Mf = aM P i anck M^ / v^2, where a is the gauge coupling and M c 
is the compactification scale. Thus one gets Mj « 10 10-12 GeV with M c « 10 8-10 GeV. 
In addition, D-brane constructions are explicit scenarios where not only intermediate 
scales arise naturally but also non-universality. 

The crucial point in these analyses is the D-brane origin of the U(l)y gauge group 
as a combination of other {7(1) 's, and its influence on the matter distribution in these 
scenarios. In particular, scenarios with the gauge group and particle content of the 
SUSY SM lead naturally to intermediate values for the string scale, in order to re- 
produce the value of gauge couplings deduced from experiments. In addition, the soft 
terms turn out to be generically non- universal. Due to these results, large cross sections 
can be obtained [0 . 



Let us consider for example a type I string scenario where the gauge group U (3) x 
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Figure 6: Scatter plot of the neutralino-proton cross section o function of 

the neutralino mass m^o in the D-brane scenario with the string scale Mj = 10 12 GeV 
discussed in the text, and for tan/3 = 10 and 15. Only the big (red and blue) dots 
fulfil b — *■ sj and g^ — 2 constraints. The red ones correspond to points with rrih > 114 
GeV whereas the blue ones correspond to points with 91 < rrih < 114 GeV. DAMA 
and CDMS current experimental limits and projected GENIUS limits are also shown. 



£7(2) x U(l), giving rise to 577(3) x SU(2) x £7(1) 3 , arises from three different types 
of D-branes, and therefore the gauge couplings are non-universal. Here £7(l)y is a 
linear combination of the three £7(1) gauge groups arising from £7(3), £7(2) and £7(1) 



within the three different D-branes. As shown in ref. |13|| , this leads to the string 
scale Mj = 10 12 GeV. Likewise, assuming that only the dilaton (S) and moduli (Tj) 
fields contribute to SUSY breaking, it was found that the soft terms are generically 
non-universal. 

Fig. H displays a scatter plot of function of the neutralino mass m^o for 

a scanning of the parameter space described in ref. fL3"l . In particular, the gravitino 
mass m 3 / 2 < 300 GeV is taken, since larger values will always produce a cross section 
below DAMA limits. The cases tan/3 = 10 and 15 are shown in Fig. Although 
regions of the parameter space consistent with DAMA limits exist, the b — > 57 and 

— 2 constraints forbid most of them. The latter are shown with small (green) 
points, and they have 91 < m h < 114 GeV. For example, to understand this result 
for ^ — 2 is simple. In ref. |34| a M versus m 3 / 2 is plotted, and one can observe that 
for m 3/ / 2 < 300 GeV too large values of a M , beyond the experimental upper bound, are 
obtained. In Fig. ^ only regions with big (red and blue) dots fulfil the above mentioned 
constraints. The red ones correspond to points with rrih > 114 GeV whereas the blue 
ones correspond to points with 91 < rrih < 114 GeV. It is worth noticing that the larger 
tan f3 is, the smaller the regions allowed by the experimental constraints become. For 
example, increasing tan/3 the value of turns out to be larger and may exceed the 
experimental bounds. 
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7 Conclusions 



We have reviewed the direct detection of supersymmetric dark matter in the light of 
recent experimental results. They require a large cross section of the order of 10~ 6 pb. 
Although in the context of the mSUGRA scenario one cannot obtain this value once 
the experimental constraints are imposed, this is possible for other scenarios. This is 
the case of SUGRA models with non-universal soft scalar masses. Likewise, scenarios 
with intermediate unification scale may also produce this value if some experimental 
bounds are relaxed. 
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